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ABSTRACT 

In astrophysical situations, e.g. in the interstellar medium (ISM), neutrals 
can provide viscous damping on scales much larger than the magnetic diffusion 
scale. Through numerical simulations, we have found that the magnetic field 
can have a rich structure below the dissipation cutoff scale. This implies that 
magnetic fields in the ISM can have structures on scales much smaller than parsec 
scales. Our results show that the magnetic energy contained in a wavenumber 
band is independent of the wavenumber and magnetic structures are intermittent 
and extremely anisotropic. We discuss the relation between our results and the 
formation of the tiny-scale atomic structure (TSAS). 

Subject headings: ISM:general—ISM:structure—MHD—turbulence 


1. Introduction 

In the interstellar medium (ISM), flows are complicated and dynamic. Observations 
suggest that the ISM is in a turbulent state (Larson 1981; Myers 1983; Scalo 1984; Armstrong, 
Rickett & Spangler 1995; Stanimirovic & Lazarian 2001). 

Hydrodynamic turbulence can be described by so-called energy cascade model, in which 
energy injected at a scale L cascades down to progressively smaller and smaller scales. Ul¬ 
timately, the energy will reach the molecular dissipation scale U and the energy will be 
lost there. The scales between L and Id constitute the inertial range. In hydrodynamic 
turbulence the dissipation scale is the minimal scale for motion. If we plot kinetic energy 
spectrum E v [k ), the kinetic energy contained in a wavenumber band, we will see that a) the 
spectrum peaks at the wavenumber corresponding to the energy injection scale [ki ~ i/i); 
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b) it follows a power law (e.g. E v (k) oc k~ 5 ^ 3 in Kolmogorov theory) in the inertial range; c) 
it drops rapidly after k c i ~ 1 /Id, which depends on viscosity v. kd ~ (KL/z/) 3 / 4 ^, where V 
is the rms velocity at the energy injection scale L. 

Magnetohydrodynamic (MffD) turbulence has two energy loss scales - a viscous damping 
scale set by the viscosity v and magnetic diffusion scale set by the ohmic resistivity r/. When 
v 3> r/, the viscous damping scale is much larger than the magnetic diffusion scale. Although 
MHD turbulence is different from its hydrodynamic counterpart in many ways, the energy 
cascade model is still valid (e.g. Goldreich & Sridhar 1995). Therefore, in this case, kinetic 
energy is damped before the cascading energy reaches the magnetic diffusion scale. In this 
paper we study MHD turbulence when the mean held B 0 is at least as large as the fluctuating 
held b. This is the opposite of the dynamo regime with B 0 <C b discussed in Kulsrud & 
Anderson (1992) and tested in Maron & Cowley (2001). We shall show that the wide-spread 
assumption that magnetic structures do not exist below the viscous cutoff is wrong. 

In astrophysics, the viscosity caused by neutrals damps turbulence. In the ISM, this 
viscous cutoff occurs at ~ pc scales, which is much larger than the magnetic diffusion scale. 
To model this, we use a large physical viscosity and very small magnetic diffusivity. 

In this letter we numerically demonstrate the existence of a power-law magnetic energy 
spectrum below the viscous damping scale. We use an incompressible MHD code. 


2. Method 

Numerical Method. We have calculated the time evolution of incompressible magnetic 
turbulence subject to a random driving force per unit mass. We have adopted a pseudo- 
spectral code to solve the incompressible MHD equations in a periodic box of size 2n (see 
more in Lazarian, Vishniac & Clio 2002): 

dw 

— = -(V x v) x V + (V x B) x B + uV\ + f + VP', (1) 

— = Vx (vxB)+ f ,V 2 B, (2) 

where f is a random driving force, P' = P/p + u 2 /2, v is the velocity, and B is magnetic 
held divided by (Anp) 1 ^ 2 . In this representation, v can be viewed as the velocity measured 
in units of the rms velocity of the system and B as the Alfven speed in the same units. The 
time t is in units of the large eddy turnover time (~ L/V ) and the length in units of L, the 
scale of the energy injection. In this system of units, the viscosity v and magnetic diffusivity 
p are the inverse of the kinetic and magnetic Reynolds numbers respectively. The magnetic 
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field consists of the uniform background field and a fluctuating field: B = Bo + b. The 
Alfven speed of the background field, B, h is set to 1. We use the same numerical technique 
as in Clio & Vishniac (2000). 

We use a physical viscosity (u = 0.015) and hyper diffusion. The power of hyperdiffusion 
is set to 2 or 3, so that the magnetic dissipation term in the above equation is replaced with 
^-r/ n (V 2 ) n v, where n = 2 or 3 and rj n (=?/ 2 or r} 3 ) is adjusted in such a way that the magnetic 
dissipation occurs around k ~ N/ 3, where N is the number of grids in each spatial direction. 
This way, we can avoid the aliasing error of the pseudo-spectral method. We use the notation 
384XY-Z?oZ, where 384 refers to the number of grid points in each spatial direction; X=P 
refers to physical viscosity; Y = H2, H3 refers to hyper-diffusion (and its power); Z=1 refers 
to the strength of the external magnetic field. 

Parameter Space. We require that the Alfven frequency is similar to the eddy turnover 
rate at the energy injection scale: B 0 k^L ~ Vk± : L, where k\\^ and k±,L are the wavenumbers 
parallel and perpendicular to B 0 at the scale, respectively. We require a large viscosity 
to maximize the dynamical range of the damped, conducting regime, but we need v small 
enough to guarantee large scale turbulence. Therefore, we take v ~ 0.015 (for a Reynolds 
number R = LV/v ~ 100). When we take v = 0.015, k L ~ 2.5, and V ~ 0.6, then 
k d = R 3/4 k L rs-/ 70. In practice, the kinetic energy spectrum usually begins to drops rapidly 
one decade before k d , so the viscous cutoff scale is around k = 7. We also require v S> rj. 
Therefore, we take hyper-diffusion for magnetic held. We consider only low order hyper¬ 
diffusion. Using hyper-diffusion, we can have a very small magnetic diffusion for small k’s 
and very strong diffusion at large k’s so that the large scale magnetic held is largely unaffected 
by ohmic dissipation. 


3. Results 

Time Evolution and Spectra. About 4 to 5 time units after the start of the simulation, 
the system has reached a statistically stationary state, where the average rms velocity (~ 0.6) 
is roughly equal to the average strength of the random magnetic held. Energy spectra reach 
a statistically stationary state a bit later at ~ 6. For the run 384PH2-Bol we integrate the 
MHD equations from t — 0 to t — 12. For the run 384PH3-B 0 1 we integrate the MHD 
equations from t — 8 to t — 12, where we use a data cube from the run 384PH2-B 0 1 at t — 8 
as the input. We calculate energy spectra and other quantities at f ~ 12. 

In Figure 1, we plot energy spectra at t=12. The spectra consist of several parts. First, 
the peak of the spectra corresponds to the energy injection scale. Second, for 2 < k < 7, 
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Fig. 1.— Energy spectra E v (k) and E b (k) at t=12. Spectra of 384PH2-B 0 1 and 384PH3-B 0 1 
almost exactly coincide at small wavenumbers. Run 384PH3-Bol has more extended inertial 
range below the viscous cutoff (marked by an arrow). 

Fig. 2.— Parallel wavenumber k\\. Solid line: 384PH2-B 0 1 at t=10. Dotted line: 384PH3- 
B 0 1 at t=12. 

kinetic and magnetic spectra follow a similar slope. This part is more or less a severely 
truncated inertial range for undamped MHD turbulence. Third, the magnetic and kinetic 
spectra begin to decouple at k ~ 7. Fourth, after k ~ 20, a new damped-scale inertial range 
emerges. We checked that the tail of the magnetic fluctuations is a real physical effect and 
not due to a bottle-neck, by comparing with calculations with a real magnetic diffusivity 
and obtaining a similar effect. 

In summary, scaling for energy spectra below the viscous cutoff is as follows. The mag¬ 
netic spectrum from Run 384PH2-B 0 1 follows 

E b (k) oc k~\ (3) 

while Run 384PH3-B 0 1 shows a slightly shallower slope. The kinetic energy spectra are 
roughly 

E v (k) oc k~ A ' 5 . (4) 

Scaling of k\\. For the kind of intermittent and extremely anisotropic structures seen these 
simulations, the interpretation of our results in terms of structure functions is complicated. 
Here we concentrate on fc||(/cj_) as a measure of field line curvature. There may be other 
independent length scales parallel to B 0 characterizing the properties of the magnetic held 
structures. 
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The term B • VB describes magnetic tension and is approximately equal to Bok\\bi for an 
isolated eddy in a uniform mean field B 0 . Here k\\ oc l/l\\ and l\\ is the characteristic length 
scale parallel to B 0 , which is known to be larger than the perpendicular length scale. Clio 
& Vishniac (2000) and Clio, Lazarian & Vishniac (2002) argued that, in actual turbulence, 
eddies are aligned with the local mean field B l. We can obtain the local frame representation 
of fc||, by considering an eddy lying in the local mean field B^: B l • Vb, m Bi,k\\bi. 1 Fourier 
transform of both sides yields |B^ ■ Vb,|k ~ |6|k, where hatted variables are Fourier- 

transformed quantities. From this, we have 
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Figure 2 shows that by this measure k\\ is nearly constant: 

k\\ ~ k ld = constant, (6) 

where k\\ td is the parallel wavenumber at the viscous damping scale. 2 In the second-order 
hyperdiffusion run 384PH2-B 0 1, k\\ is almost constant up to k ~ 50. The sharp rise after 
k ~ 50 is due to magnetic dissipation and the exponential suppression of the magnetic field 
energy. The third order hyperdiffusion run (the dotted line) shows more extended range of 
constant k\\ since magnetic dissipation occurs at larger /c’s. 

We checked whether or not equation (5) gives a reasonable result for ordinary (i.e. 
not viscously damped) MHD turbulence. Clio & Vishniac (2000) showed numerically that 
k\\ oc k 2 / 3 using structure functions when the turbulence is threaded by a strong external 
magnetic field (see Goldreich & Sridhar 1995 for a theoretical derivation). When we applied 
the method in equation (5) to this case we recovered the relation k\\ oc k 2 / 3 , confirming 
our intuitive notion that equation (5) can provide us with a coordinate-independent way of 
calculating k\\ as a function of k±. 

Intermittency. In Figure 3, we show strength of small scale magnetic fields in a plane 
perpendicular to the mean field Bo- Darker tones represent stronger magnetic fields. We also 
plot small scale magnetic vectors projected on the plane. We obtain the small scale magnetic 


Where can be many ways to define the local mean field B^. In this paper, we obtain Bl for an eddy of 
(perpendicular) size l oc 1/k by eliminating the Fourier modes whose perpendicular wavenumber is greater 
than k/2 and b, by eliminating the Fourier modes whose perpendicular wavenumber is less than k/2. 

2 Constancy of k || means an extreme form of anisotropy: eddies have a constant parallel size while they 
can have extremely small perpendicular size. 
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field by eliminating Fourier modes with k < 20. Only part of the plane is shown here. We 
see that the magnetic structures are very intermittent. 3 We also see that the typical radius 
of curvature of field lines in the plane is much larger than the typical perpendicular scale 
for field reversal. The typical radius of curvature of field lines corresponds to the viscous 
damping scale, indicating that stretched structures are results of the shearing motions at 
the viscous damping scale. This is consistent with Figure 2, which shows that the parallel 
wavelength at the viscous damping scale is also reflected in the curvature of the small scale 
field lines. There is no preferred direction for these elongated structures. A similar plot for 
ordinary (not viscously damped) MHD turbulence shows much less intermittent structures. 

What is the Scaling? In this section, we give a simple set of scaling relations for the 
properties of this regime (e.g. 20 < k < 50 for Run 384PH3-B 0 1). Later, in Lazarian, 
Vishniac & Clio (2002, henceforth LVC02), we will discuss the physical basis for these scaling 
relations in detail. 

We first note that the magnetic spectrum is roughly proportional to k~ l . When the 
shearing by the motion at scale l ( i dominates that of the scale /, the damping-scale shearing 
timescale ( ld/v<i ) becomes the characteristic timescale at the scale l. Note that the charac¬ 
teristic timescale is scale-independent. Assuming a local cascade of energy in phase space 


we have 


bf/(ld/vd) ~ constant, 

(7) 

so that bi ~ constant, which, in turn, implies 


E b (k) oc k~\ 

(8) 


This is identical to the viscous-convective range of a passive scalar in hydrodynamic turbu¬ 
lence (see, for example, Lesieur 1990). 

Figure 2 suggests that k\\ ^constant below the cutoff. 

Figure 3 shows that the magnetic structures are very intermittent below the viscous 
cutoff. Assuming that the filaments are characterized by a rough balance between viscous 
drag and magnetic tension, calculations in LVC02 predict that 

E v (k ) oc k~ A (9) 

which is almost exactly true in our simulations. 


3 This is a real physical effect. When we use a real magnetic diffusivity, we also obtain highly intermittent 
structures similar to those in Fig. 3. 



4. Implications 


The intermittent small scale structures that we predict here should have important im¬ 
plications for transport processes (heat, cosmic rays, etc.) in partially ionized plasmas. We 
also speculate that they might have some relation to the tiny-scale atomic structures (TSAS). 
Heiles (1997) introduced the term TSAS for the mysterious H I absorbing structures on the 
scale from thousands to tens of AU, discovered by Deiter, Welch & Romney (1976). Analogs 
are observed in Nal and Call (Meyer & Blades 1996; Faison & Goss 2001; Andrews, Meyer 
& Lauroesch 2001) and in molecular gas (Marscher, Moore & Bania 1993). Recently Desh- 
pande, Dwarakanath & Goss (2000) analyzed channel maps of opacity fluctuations toward 
Cas A and Cygnus A. They found that the amplitudes of density fluctuations at scales less 
than 0.1 pc are far larger than expected from extrapolation from larger scales, possibly 
explaining TSAS. This study, however, cannot answer what confines those presumably over¬ 
pressured (but very quiescent!) blobs of gas. Deshpande (2000) related those structures to 
the shallow spectrum of interstellar turbulence. 

Figure 1 indicates that while velocity decreases rapidly, but not exponentially, below the 
viscous damping scale, the magnetic field fluctuations persist, thereby providing nonthermal 
pressure. Magnetic structures perpendicular to the mean magnetic field are compensated 
by pressure gradients from density fluctuations reminiscent of the Dishpande et al. (2000) 
observations. 

Our calculations are applicable on scales from the viscous damping scale (determined 
by equating the energy transfer rate with the viscous damping rate; ~ 0.1 pc in the Warm 
Neutral Medium with n = 0.4 cm -3 , T— 6000 K) to the ion-neutral decoupling scale (the 
scale at which viscous drag on ions becomes comparable to the neutral drag; <C 0.1 pc). 
Below the viscous scale the fluctuations of magnetic field obey the damped regime shown in 
Figure 1 and produce density fluctuations. If gas is close to be thermally unstable (Vazquez- 
Semadeni, Gazol & Scalo 2000), moderate fluctuations of pressure cause substantial density 
excursions. For typical Cold Neutral Medium gas, the scale of neutral-ion decoupling de¬ 
creases to ~ 70AU, and is less for denser gas. TSAS may be created by strongly nonlinear 
MHD turbulence! 

The issue of an adequate description of turbulence in the newly found regime of per¬ 
sistent turbulence goes much beyond the TSAS formation. The accumulation of energy at 
small scales is important for problems of cosmic ray transport, magnetic reconnection, and 
formation of clumps in molecular clouds. We expect to make significant progress in those 
fundamentally important directions as soon as we describe properly the new regime of MHD 
turbulence. 



5. Conclusion 


We have considered MHD turbulence with v 3> r/, which implies that the viscous cutoff 
occurs at a scale much larger than the magnetic diffusion scale. It has been believed that the 
damping of the fluid motion is accompanied by the suppression of magnetic structures below 
the viscous cutoff scale. To the contrary, we have found that the magnetic held can have a 
rich structure below the viscous cutoff scale. Judging from our (limited) numerical results, 
we conclude that magnetic held perturbations have a similar power distribution as a passive 
scalar, despite the obvious importance of magnetic forces. In particular, the spectrum follow 
a k~ x law. The parallel wavenumber, which is an indicator of the degree of anisotropy, is 
almost constant. Consequently, the magnetic held has an extreme form of anisotropy. In 
summary, we have found that 

1. Eb(k) oc /c _1 , 

2. E v (k ) oc /c -4 , 

3. fc|| ~ constant. 

We discussed the possibility that this small scale magnetic held is the cause of the tiny-scale 
atomic structure. The small scale magnetic structure will affect many astrophysical processes 
(e.g. cosmic-ray transport, reconnection, etc) that depend on the statistical properties of 
MHD turbulence. 
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